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I. INTRODUCTION 


The Naval Undersea Warfare Engineering Station (NUWES) conducts testing of 
torpedoes on various underwater ranges in the Pacific Northwest (Washington State and 
Canada). On occasion these torpedoes malfunction and are lost on or under the ocean 
floor and must be recovered. The recovery operation involves locating the torpedo, if 
neccessary digging it out from under the ocean bottom, and attaching a cable to it in 
order to pull it out. This is all done by remote control by an operator on a surface ship 
using video and sonar images to guide him. The sonar signal is used for general position 
location and then video images are used to conduct the remainder of the operation. 

There are several problems which occur during the acquisition of these images. 
Probably the most critical of these are caused by the artificial lighting required at depths 
of up to 1500 feet. Uneven lighting can produce severe shadows and bright blooming 
on the same image frame. Bottom texture (important for finding the entry point when 
the torpedo has buried itself) often cannot be distinguished. Silt and debris thrown up 
by the digging device can totally obscure vision for up to several hours at a time. 

The tools chosen to help solve these problems were an IBM AT Personal 
Computer! with a PCVISIONplus FRAMEGRABBER2? board installed and the ITEX 
PCplus3 library of image processing subroutines. This svstem allows one to build and 
run software tailored to his her particular image processing needs. The svstem 1s rela- 
tively inexpensive (when compared to more sophisticated image processing work- 
stations), portable, and easily installed in an office or on board a recovery vessel. 

The purpose of this thesis was to develop an image processing package, using the 
tools described above, to provide contrast and other enhancement of video images. Be- 
cause the images to be processed are used as feedback in operating a remote control 
device, a real-time or near real-time approach to the problem was neccesitated. Ina 
previous thesis, [Ref. 1] Roberto Ventura suggested the use of lookup table modification 
as a primary solution. There were also several other non-lookup table algorithms which 


appeared to have potential usefulness. Once developed, the package was to be tested 


1 Trademark of IBM 
2 Trademark of Imaging Technology Inc. 


3 Trademark of Imaging Technology Inc. 


at NUWES and delivered to them for their use. Inherent in this requirement was a 
“user-friendly” presentation at the computer monitor and appropriate documentation for 
its operation. The text of this thesis shows an analvsis of the problem, describes, in 
general, the tools used, and outlines the development of the final package. Further, it 
describes the results of the testing conducted on a recovery vessel at NUWES. The ap- 
pendices include instructions for the use of the package and the C language code for the 


program. 


to 


II. ANALYSIS OF THE PROBLEM 


After observing a recovery operation, receiving comments from operators, and ex- 
amining video tape of typical underwater conditions, three viewing regimes were identi- 
fied. These were viewing through turbid water (due to disturbance of bottom silt), clear 
water viewing of the bottom while searching, and clear water viewing of shiny objects 
which may cause blooming. Each regime was analyzed to determine what problems, if 
any, could potentially be solved. The calculation and display of histograms, a primary 
tool in the analysis, will be discussed in a later chapter. 

The most challenging regime is that of viewing objects through turbid water. Often 
the densitv of particles in the water is so great that all incident light is scattered or re- 
flected before it can reach the area or object to be observed. In this case the video sys- 
tem is of no use and some other portion of the electromagnetic spectrum should be 
considered for image acquisition. However, if even a little light is reflected from the 
object of interest back to the camera there 1s potential for enhancement of the resulting 
image. The histogram of a typical turbid water image, as shown in Figure 1, provides 
information which could lead to the method used to acheive this improvement. The 
tight distribution of pixel values around the center range of intensity levels readily lends 
itself to contrast enhancement using look-up table modification techniques. 

When searching for a torpedo that has buried itself under the bottom, the general 
location is first determined using sonar. Once in the vicinity, however, it is verv useful 
if the entry point can be found. This is often difficult since the mud closes in over the 
top of the torpedo. It would be helpful if some method could be found to better identify 
marks on the bottom which could be entry points. Figure 2 shows the histogram of a 
tvpical image of the bottom. Again, most pixel values are concentrated around the 
center range of intensity levels and this provides excellent opportunity for enhancement 
by look-up table modification. Also, an enlarged view of a portion of the image could 
be potentially helpful in object identification. This could aid in identifying a partially 
uncovered object during the digging operation. 

Blooming is a phenomona caused by the specular reflection of bnght lights from 
shiny objects. It often occurs after the torpedo is found while attempting to attach a 
cable to it. The result is a “white out” over a significant portion of the image. The 


current Operational solution to this problem is to dim the entire image on the monitor 





Figure 1. Histogram of a Turbid Water Image 





Figure 2. Histogram of a Bottom Image 


using the brightness control. This reduces the detail available from the darker portions 
of the image. The histogram of an image with blooming, shown in Figure 3, also sug- 


gests that a look-up table modification could be performed to reduce the brightness of 


higher intensity pixels without dimming the lower values. 





Figure 3. Histogram of Image With Blooming. 


II. DESCRIPTION OF THE IMAGE PROCESSING SYSTEM 


A. GENERAL 

This chapter describes the tools to be used in the attempt to improve the quality of 
the underwater images. The image processing system, shown in Figure 4, includes 
hardware in the form of a _ personal computer with the ITEX PCplus 
FRAMEGRABBER installed, and a library of image processing software called ITEX 
PCplus. 


B. HARDWARE 

The framegrabber (enclosed within the dashed line in Figure 4) is the heart of the 
image processing system. It is installed in a host computer for control of its components 
and for disk access. Compatibility with the host is ensured by placement of jumpers in 
the proper locations on the board. This also sets the memory base and register base 
addresses. The framegrabber receives video input from a source (camera or video tape 
plaver) and sends video output to a monitor. At the input the image is digitized by an 
analog-to-digital converter. It is then sent through the input lookup table (LUT). There 
are eight hardware lookup tables available on the input side. One 1s selected for use at 
any given time and can be loaded with the desired values from the computer. Values 
mav also be read from the selected lookup table and stored to disk. After passing 
through the input lookup table the image is placed into one of two memories, A or B, 
whichever is selected from the computer. The contents of these memories are accessible 
from the computer for analvsis, processing, or storage to disk. Once an image is in a 
selected memory, it may be sent for display. Enroute it passes through one of eight 
output lookup tables which operate in the same manner as the input lookup tables. 
Finally it is converted to an analog image to be displayed on the video monitor. In the 
case of a continuous video image from the source this entire process may occur every 
one-thirtieth of a second, presenting a continuous image to be displayed at the output. 
Single frames from a source or from the computer disk mav also be placed into memory 


and presented at the output for display. [Ref. 2] 


C. “SOFTWARE 
The framegrabber hardware is controlled using a library of functions called ITEX 


PCplus. These functions may be incorporated in a program using the “C” language. 
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Figure 4. The Image Processing System 





They perform various tasks including computer setup, image acquisition, image Storage 
to disk, lookup table operations, and single frame image processing. Along with basic 
functions to place the framegrabber into operation, several others were identified as 
having potential for providing real time or near real time contrast enhancement. These 
building blocks are described here. 

The setup functions configure the computer to be compatible with the framegrabber 
configuration, select memory size options, initialize the framegrabber, and place it into 
operation. They also allow selection of memory for acquisition, display, storage, and 
processing. A brief description of these functions follows: 


sethdw Defines the hardware settings including register base address, memory base 
address, and memory option. Register base and memory base address ar- 
guments tell the computer how the framegrabber is configured. The mem- 
ory option defines either dua] 512x512 frame memories or a single 1024x512 
memory. 


setdim Defines the memory size. Sets memory’s horizontal and vertical size and 
number of bits used to store each pixel. 


initialize Initializes all hardware. Sets all hardware registers and lookup tables to 
standard values. 


foon Turns on the framegrabber memory. 


select_mem Selects frame memory. Allows selection of memory for loading from disk, 
storage to disk, or processing. May be used to access one memory while 
another memory is being displayed. 


display_mem Selects frame memory display. Allows selection of memory for acquisi- 
tion from source and/or display to screen. 


feoff Turns off the framegrabber memory. [Ref. 3] 


The image acquisition functions allow input of image data from the source. They 
allow continuous passage of image frames from source to video monitor, single frame 


storage and display, and memory clear operations. They are described as follows: 


snap Acquires a single image frame from the source, stores it in frame memory, 
and displays it on the monitor. Allows single frame processing or storage 
to disk of an image from the source. 


grab Acquires and displays image frames continuously up to a specified number 
of frames or indefinitely until new a command 1s given. Allows real time 


image processing since image frames pass continuously through lookup ta- 
bles. 


sclear Clears entire screen to a value. Erases frame memory and screen by placing 
a specified intensity level into each memory location. 


saveim Saves a specified portion of an image to a file on disk. The file name 1s 
passed as an argument to the function. 


readim Reads an image file from disk. The file name 1s passed as an argument to 
the function. [Ref. 3] 

The lookup table functions provide the capability to select a specified lookup table 

for use, read lookup table values from disk, write lookup table values to disk, and modif¥ 


lookup tables in various wavs. These functions are described below: 


setlut Selects which of eight lookup tables on the input or output bank will be 
active and available for storage of values. 

walut Writes an array of values to a selected lookup table. 

ralut Reads selected lookup table values into an array. 


save_lut Saves the values in a selected lookup table to disk. 


read_lut Reads lookup table values from disk and places them into the lookup table 
from which thev were stored. 


linlut Places linear lookup table values into selected lookup table. 


invlut Places inverse linear lookup table values into selected lookup table. 


loglut Places logarithmic lookup table values into selected lookup table. 


abslut Places values into selected lookup table which invert all inputs greater than 
127. 
histeq Performs histogram equalization on a specified portion of an image and 


places computed values into selected lookup table. [Ref. 3] 


Several non-real time image processing functions were also considered. These 
functions work on images fixed in frame memory and perform filtering, pixel replication, 
averaging, and histogram operations. They are described here: 


sobel Performs a Sobel edge detector operation on a specified portion of an im- 
age. Displays prominent edges of the image. 


sharpen Sharpens a specified portion of an image by selectively amplifying edges. 


lopass Performs a low-pass filtering operation on a specified portion of an image 
and displays result. 


zoom Enlarges by a factor of two the portion of the image in the upper left 
quarter of the displav. 


average Acquires the selected portion of a specified number of live image frames, 
stores them on disk, computes the average, and displays results. 


histogram Computes a histogram of a specified portion of an image and places values 
into an arrav. Allows selection of spatial sampling rate in vertical and 
horizontal directions and size of storage bins (array). 


pan Shifts the entire image horizontally by a specified amount. 
scroll Shifts the entire image Vertically by a specified amount. 


rectangle Draws a rectangle with specified size and location on the screen and in 
frame memorv. May be superimposed on a fixed image. [Ref. 3] 


IV. DEVELOPMENT OF THE MENU PROGRAM 


A. GENERAL 

The concept of the “Menu” program was to place the various image processing 
options and the analytical tools into a format which allows option and parameter 
changes to be easily made and evaluated. Further, the menu format is used in the final 
image processing package to be delivered to NUWES. 

This chapter describes how the various ITEX PCplus functions were incorporated 
using Microsoft4 C 5.0 [Ref. 4] to produce the “Menu” program. It outlines how the 
setup functions were used to bring the board on line and perform some basic operations. 
Then it describes how the histogram and lookup table read’write functions were used to 
build the analvtical tools. Finally, it shows how various lookup table modification and 
other image processing functions Were written into the program and tested for applica- 


bility to the stated problems. 


B. BASIC OPERATIONS 
1. Setup 
Four functions are required to define and initialize the hardware and place the 
board into operation. The sethdw function is used with arguments which define the 
memory base address as DOOOO., the register base address as 100, and select the dual 
memory model. It 1s also used in a menu option which allows the operator to interac- 
tivelv enter different addresses if the jumpers on the framegrabber board are configured 
differently. The setdim function is used to set the frame memory size to 512 x 512 pixel 
locations with § bits for each pixel value. The initialize function 1s included in the pro- 
gram to set all registers and lookup tables to standard values and fgon 1s included to turn 
the framegrabber memory on. 
2. Image Acquisition 
Several menu options were devoted to image aquisition functions. Two options 
use the display_mem function to allow the operator to choose which frame memory to 
display. Iwo other options use the select_mem function to allow selection of a frame 
memory for processing or storage to disk. Combinations of these four choices permit 


one image to be displaved while another is being processed. The snap function was used 


4 Trademark of Microsoft Corporation 
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in a menu option to allow the acquisition of a single image frame into a selected mem- 
ory. Another option selects continuous acquisition and display of image frames using 
the grab function. The clear screen menu options use the sclear function to erase the 
screen and frame memory to either black or white. Finally, access to images stored on 


disk is accomplished with menu options using the saveim and readim functions. 


C. BUILDING THE ANALYSIS TOOLS 
1. Computation and Display of the Histogram 

The histogram was the most important tool used in analysing the characteristics 
of an image. The histogram function provides the computed histogram values in an ar- 
ray, but methods needed to be devised to specify the portion of the image to include in 
the calculation, to display the histogram graphically, and to compare images and their 
histograms before and after processing. 

To specifv the portion of the image to be used in the histogram calculation, a 
rectangle whose size and position could be changed from the keyboard was superim- 
posed on a fixed image displayed on the monitor using the rectangle function. This se- 
quence (the snap and rectangle functions) was placed inside a loop to give the 
appearance of continuous acquisition and displav. The size and location parameters of 
the rectangle were then used as the area arguments for the histogram function. A default 
option to use the entire screen was also programmed so that the rectangle routine could 
be bypassed if desired. 

A method was also needed to display the histogram values graphically on the 
computer monitor at the same time the image was being displayed on the video monitor. 
This was accomplished in a routine which used Microsoft C 5.0 linedrawing and other 
graphics functions to draw vertical lines from a horizontal axis up to scaled histogram 
values. The result was a 256 line wide (one line for each intensity level) bar graph 
showing the relative density of each intensity level of the image. Examples of the 
histogram displav are shown in Chapter 2 (Figures 1, 2, and 3). 

To make a comparison of an image and its histogram before and after process- 
ing requires the use of both frame memories. An image frame is acquired into one 
memory using the snap function with a linear input lookup table. The other memory 
is selected (the display_mem function) and new values are loaded into the input lookup 
table using the walut function. Then another image frame is acquired. The time differ- 
ence between the two frames is so short that the two images are virtually the same. The 


histogram values for each frame are then computed with the histogram function and 


stored in two arrays. The histogram and frame memory to be displaved can then be 
selected from the keyboard, presenting processed and unprocessed versions of the image 
and histogram to be examined. 

2. Manual Modification and Display of Lookup Tables 

As stated earlier, the lookup table is the primary means by which the contrast 
of the image frames can be enhanced in real time. In order to predict and evaluate the 
effect of a particular set of lookup table values on an image a means of displaying those 
values graphically was needed. Further, it was desired to be able to change the lookup 
table graph interactively from the computer keyboard and evaluate the effect of these 
changes on the image frames. 

To plot the lookup table the ralut function 1s used to read the values from the 
selected lookup table into an array. Microsoft C 5.0 linedrawing functions with line 
segment endpoints defined by arrav values plot a continuous line which represents the 
lookup table. Figure 5 shows a typical lookup table display. 

Two methods were programmed to interactively modify lookup tables. The first 
method reverses the process used to plot the lookup table described above. Line seg- 
ment end points are entered from the computer keyboard. The line is then drawn on the 
computer monitor. The equation for that line 1s computed and values are derived to be 
placed in the appropriate positions in the lookup table array. The walut function then 
writes the array to the lookup table. The effect on the image is immediately evident and 
may be evaluated. The second modification routine interactively varies the slope of the 
lookup table. The point where the line intercepts the top of the graph 1s shifted to the 
left or right bv pressing the arrow keys on the computer kevboard. The remaining points 
on the line to the bottom of the graph are then shifted in the same direction, with the 
amount of shift varving linearly to zero from the initial shift. The shape of the line re- 
mains very close to the original, but the overall slope is changed, as shown in Figures 6 
and 7. As the point by point changes are made, the lookup table array is continually 
updated and the values written to the selected lookup table. The effect of each change 
on the incoming image frames occurs immediately. The slope may also be varied in a 
similar manner bv shifting the bottom portion of the lookup table line. 

At this point the save_lut and read_lut functions were included as menu options 
to allow future use of particularly effective sets of lookup table derived during exper- 


imentation. 





Figure 5. Typical Lookup Table Display 
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Figure 6. Lookup Table Before Slope Modification 


OUTPUT 





Figure 7. Lookup Table After Slope Modification 
D. APPLICATION OF LOOKUP TABLE OPERATIONS 
1. Lookup Table Stretch Routines 

As shown in Chapter II of this thesis (Figures 1 and 2), the histograms for 
turbid water images and bottom images have high concentrations of pixels in the center 
range of intensity levels and few or zero in the outer ranges. In these situations contrast 
enhancement may be achieved by using lookup table stretch methods to expand the 
histogram to cover the full range of available intensity levels. [Ref. 5] The following 
experiment was performed. The histogram (shown in Figure 8) of a tvpical image from 
of a sequence of bottom images was evaluated visually for the high and low intensitv 
values of the high pixel concentration area. Then a linear lookup table was modified 
using the slope varying routine described in this chapter so that the top intercept 
matched the high end of intensity level concentration and the bottom intercept matched 
the low end of intensity level concentration as shown in Figure 9. The histogram routine 
was then apphed to compare the enhanced image with the unenhanced image. The 
modified histogram is shown in Figure 10. The result was an image whose contrast Was 
significantly enhanced. The problem then became one of automating the histogram 
evaluation process and making it responsive to changes as they occurred in the contin- 


uously acquired images. 





Figure 8. Histogram Before Lookup [able Stretch 
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Figure 9. The Lookup Table Stretch 





Figure 10. Histogram After Lookup Table Stretch 


Histogram equalization provides the answer to the histogram evaluation portion 


of the problem. This process is defined bv a transformation function of the form 


Vix) = | “pls 


where p/s) 1s the histogram distribution function, s is a dummy variable of integration, 
x 1s a normalized input gray level, and V is the transformed output gray level. This 
transformation function, when plotted for the full gray scale, gives the shape of the de- 
sired lookup table. [Ref. 6] The histeq function uses this concept to find the output 
value for each input to the lookup table by computing the normalized cumulative area 
up to that input intensity level on the histogram. It then loads those output values into 
the specified lookup table. In effect, the function steepens the slope of the lookup table 
plot for high concentration portions of the histogram and shallows the slope for low 
concentration portions, a localized lookup table stretch. Figures 11 and 12 show the 
lookup table after equalization and the resulting histogram. Figures 13 and 14 show an 
image before and after equalization. Although the overall image appears darker, details 


of the bottom are more readily apparent. 
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Figure 11. Lookup Table After Equalization 





Figure 12. Histogram After Equalization 
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The problem of changing image frame content and hence changing histogram 
was solved by programming two alternatives for update (performance of a new 
histogram equalization). The first is a manual update performed from the computer 
keyboard bv the operator as he notices the image content changing on the video mon- 
itor. The second is an automatic update where the computer routinely grabs a frame 
and performs a new equalization at a specified time interval. With either alternative the 
update occurs very quickly with less than a one second break in the continuous image 
acquisition as the new lookup table 1s loaded. 

An option for equalization using the intensity levels of a specified portion of the 
image was also included. The user designates the portion of the image using the same 
rectangle described in the histogram portion of this chapter. The rectangle parameters 
are used to define the area used for the equalization; this results in maximum contrast 
for that area. 

2. Other Lookup Table Routines 

Chapter IJ of this thesis also describes the problems caused by the lighting re- 
quired in the underwater environment. The histogram (shown in Figure 3) of an image 
in which blooming was occurring suggested that a modification to the portion of the 
lookup table containing the high intensity levels was in order. The lookup table drawing 
routine described in this chapter was used to experiment with different lookup table 
shapes which reduced the intensity of high value pixels. An example of one of these is 
shown in Figure 15. Also, the invlut and abslut functions were tried. These lookup table 
shapes are shown in Figures 16 and 17. The procedures most effective in reducing the 
blooming without affecting the remainder of the image used the inverse and absolute 
value lookup tables. Both were included as menu options. 

Finally, the loglut function was tested on images which appeared dark due to 
shadows. The logarithmic values had a very definite brightening effect on the images, 
similar to that of histogram equalization on the same sequence. The loglut function was 
included as a menu option as a simple, one keystroke alternative to equalization for dark 


images. The logarithmic lookup table is shown in Figure 18. 


E. APPLICATION OF OTHER OPERATIONS 

Several non-real time image processing functions were found to be effective in en- 
hancing turbid water and bottom images and were included as menu options. These 
functions all operate on a single image frame fixed in memory. The time required to 


perform these operations varies from four to thirty seconds depending on the operation 
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Figure 15. Experimental Lookup Table 
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Fieure 16. Inverse Lookup Table 
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Figure 17. Absolute Value Lookup Table 
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Figure 18. Logarithmic Lookup Table 
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and computer speed. Analvsis of the effectiveness of these options was limited to visual 
evaluation of the actual image. 

The zoom function was included as a menu option to assist in object identification. 
In order to allow selection of the portion of the screen to be enlarged, a rectangle routine 
similar to that used in the histogram and equalization options was programmed. In this 
case, however, the rectangle location parameters are sent to the pan and scroll functions 
to shift the image in memory to the upper left quadrant, the portion of memory used in 
the zoom function. Once the image is enlarged, it may be processed by other options 
that operate on a single frame fixed image. 

The sharpen function was included as an option to provide the capability to amplify 
the edges of a fixed image. It is potentially helpful in object identification or in finding 
an entry point on the bottom. 

The lopass function provides the capability to filter out high frequecy noise from a 
fixed image. It is included as a menu option as a potential tool for turbid water viewing. 

The sobel function was included as a menu option to provide an edge detection ca- 
pability. It is potentially useful in object identification. 

Image averaging was included as a menu option using the average function. When 
this option is called, a specified number of image frames are acquired and stored to disk. 
They are then averaged, pixel bv pixel, and the result 1s displaved on the video monitor. 
This option has potential usefulness when operating in turbid water with the camera 
mounted on a stable platform. It allows the selection of the number of image frames to 
be averaged up to disk memory capacity. Itis the slowest of the processing options due 
to the high number of read write to disk operations and large number of calculations 


required. 
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V. ON SITE TESTING 


A. GENERAL 

An important aspect of this project was to provide the NUWES recovery section 
with a package that was easy to use and helpful in their specific application. For this 
reason a dialogue was established and maintained with the NUWES engineers and re- 
covery section operators for the duration of the development process. The author made 
three visits to NUWES. The purpose of the first visit was to observe a recovery opera- 
tion first hand and define the key elements of the problem. The next visit was conducted 
midway into the project. Its purpose was to demonstrate a prototype “Menu” program 
and obtain feedback about the usefulness of the various image processing options and 
the format of the on-screen presentation. The final visit was an on-site test of the pro- 
gram, conducted on a torpedo recovery vessel with cameras and equipment on the bot- 
tom under actual recovery conditions. The aft section of a torpedo was used as the test 
object to be viewed. The results of that test and the feedback obtained from NUWES 
operators and engineers are the subject of this chapter. 


B. SPECIFIC APPLICATIONS 

The three specific applications initiallv identified for processing were turbid water 
viewing, bottom viewing while searching, and blooming. Additionally, during the course 
of the development, enhancement of shadowed images and object identification became 
areas of interest. 

The operations expected to be most effective during turbid water viewing were 
histogram equalization, the low pass filter, and image averaging. This portion of the test 
provided the most difficult environment. The histogram equalization did little to im- 
prove the the overall quality of viewing. The periods when the object was visible 
through the silt were so short that by the time a manual update was performed the ob- 
ject had disappeared. Similarly, in the automatic update mode the randomness of the 
short appearances of the object defeated the periodic update of the lookup table. The 
situation was further complicated by a horizontal dark grav synchronization bar moving 
through the image every two to three seconds. This reduced the level of effectiveness 
of the equalization since the bar’s intensity values were included in the computation. 
In fact, the equalization process enhanced the bar. The svnchronization bar had not 


been visible during the previous on board visit or on videotape and so may have been a 


problem peculiar to the particular equipment used on this test. The low pass filter also 
proved to be of limited value in this environment. Again, the periods when the object 
was visible at all were so short that it was practically impossible to freeze an image to 
be processed. 

Image averaging met with the highest degree of success during the turbid water 
portion of the test. With the recovery device (camera platform) sitting firmly on the 
bottom, the image average routine was able to capture those fleeting moments when the 
object was Visible along with the time when it was not. The averaging computation then, 
to some extent, removed the turbidity from the image. The result was a non-real time 
image in which the object was recognizable. The time elapsed from the beginning of the 
process, averaging six frames, Was approximately one minute. This was attributed to a 
relatively low computer speed (6 megahertz) and the read/write operations to the hard 
disk.5 

The operations included in the program for enhancement of images during bottom 
viewing Were histogram equalization and image sharpening. During the test a camera 
was pointed at the bottom and these operations were performed. The histogram equal- 
ization routine produced good results in increasing the contrast, making bottom features 
more easily identifiable. The synchronization bar again reduced effectiveness, but the 
value of the equalization was proven for this application. The sharpen routine was also 
helpful in the enhancement of bottom images, although this routine does not operate in 
real time. It was concluded that both of these routines would be of considerable value 
in identifving torpedo entry marks on the bottom. 

The inverse and absolute value lookup table operations were included in the “Menu” 
program to help solve the problem of blooming. During the test an attempt was made 
to induce blooming by directing the lights towards the test object. However, due to the 
lighting angle and relatively dull surface of the object, severe blooming as seen on vide- 
otape did not occur. The inverse and absolute value operations performed as expected 
on the minimal blooming induced, and operator feedback indicated that these operations 
would have some utility during an actual recovery. Interest was also expressed in the 
Capability to experiment with different lookup table shapes using the lookup table 


drawing routine and save those which helped solve the blooming problem. 


> An attempt to speed up the umage averaging operation 1s descnbed in part C of this chapter. 


The zoom operation was also demonstrated and although there was no specific ap- 
plication for its use during the test, it was agreed that it would be helpful in certain in- 


Stances during recovery operations. 


C. USER FEEDBACK 

Upon completion of the test, comments were solicited from the recovery section 
operators and engineers concerning ways to improve the “Menu” program. The sug- 
gestion was made that the program’s menu presentation could be revised to allow the 
use of a mouse to make selections and move rectangles and pointers on the video mon- 
itor. A high degree of interest was also expressed in improving the image averaging 
routine. Methods suggested included compensating for motion of the camera platform, 
using weighted averaging techniques, and speeding up the process. Once ashore, an at- 
tempt to provide more speed was made by installing the framegrabber board on a 20 
megahertz Compaq 386 computer and copying the “Menu” program to a virtual disk. 
The combination of higher computer speed and use of RAM for temporarv image stor- 
age produced a significant reduction in the time required to perform the operation (ten 
frames were averaged and the result was displayed in ten seconds). 

Another area of interest was the potential to use the color capabilitv available by 
using three framegrabber boards. Suggestions ranged from coloring the rectangle 
superimposed on image frames during some operations to creating pseudocolor images. 

In general, there was enthusiasm expressed for the use and continued development 


of the capability provided by this image processing system. 
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VI. CONCLUSIONS 


The problem of enhancement of underwater images during torpedo recovery oper- 
ations was analysed and broken down into processing for three specific viewing envi- 
ronments. These were viewing through turbid water, viewing of bottom features, and 
viewing when blooming occurs. The tools available to solve these problems were then 
studied. These included an IBM AT Personal Computer with the PCVISIONplus 
FRAMEGRABBER installed and a library of image processing functions called ITEX 
PCplus. The “Menu” program was then developed by incorporating the image process- 
ing functions into an interactive program with a menu format. The program provides 
convenient and easv access to the image processing functions provided by ITEX PCplus. 
It also provides capability for automatic or manual equalization, a segmentation routine 
for improved contrast on a selected portion of an image, graphical input of lookup ta- 
bles, modification of existing lookup tables, and storage of lookup tables to disk. The 
histogram and lookup table display routines were built and included as analysis tools. 
These provide on screen display of lookup tables and histograms and the capability to 
display enhanced and unenhanced versions of an image and their respective histograms 
for comparison purposes. Some non-real time functions included were image averaging, 
a Sobel edge detector, a low pass filter, a sharpening routine, and a zoom feature which 
allows interactive selection of the portion of an image to be enlarged. Documentation 
for the “Menu” program was written in the form of a User’s Manual (see Appendix A) 
which includes framegrabber installation instructions, a tutorial on the use of the 
“Menu” program, and a quick reference list. 

A test of the final package was conducted aboard a recovery vessel at NUWES with 
cameras and equipment on the bottom under actual recovery conditions. The program 
was found to be particularly effective when viewing bottom features and during bloom- 
ing. The equalization procedures, lookup table drawing capability, absolute value 
lookup table modification, image averaging, and the zoom function were all determined 
to be potentially valuable tools for actual recovery operations. Appreciation was also 
expressed for the documentation provided with the program. 

More development was called for in the turbid water viewing environment, specif- 
ically with regard to the improvement of the image averaging routine. The current 


function averages a finite number of image frames chosen by the user, each weighted 
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equally. All the frames are stored on the hard disk and then read back into RAM for 
averaging. The speed of the process 1s limited by the clock speed of the computer and 
the read; writes to and from disk. A stable camera platform is essential or the resulting 
image becomes extremely blurred. Improvements could include continuous weighted 
frame averaging, compensating for camera motion by cross-correlating image frames, 
and use of a faster computer with a virtual disk for increased speed. 

The “Menu” program with documentation was delivered to NUWES for operational 


use. 


APPENDIX A. USER’S MANUAL 


A. GENERAL 

The “Menu” software package is a menu-driven program which provides real time 
continuous video image processing using lookup table modification techniques. It also 
provides limited non-real time single frame enhancement capability. Control of the 
routines in “Menu” varies from manual (requiring operator input to make an adjust- 
ment) to automatic or “hands-off” operation. The processing routines were chosen to 
provide maximum effectiveness in situations where the images’ pixel values are limited 
to a relatively small region of the total gray scale range. Typically this would include 
images of objects which have actual small variations in shading or those whose contrast 


has been reduced by obscuration or dim lighting. 


B. EQUIPMENT REQUIRED 
The following equipment is required to use the “Menu” program: 


¢ IBM Personal Computer AT, XT, or PC6 or 100% compatibles with Extended 
Graphics Adaptor card and monitor. 


e Imaging Technology PCVISIONplus Frame Grabber? card with standard cable. 
(NOTE: The older model PCVISION Frame Grabber will not work.) 


e Analog video monitor (such as SONY Trinitron8’ Model 1271Q). 


e Video source (camera or video tape plaver). 


C. SIMPLE SETUP 
The PCVISIONplus Frame Grabber User’s Manual provides all the information 
required to place the Frame Grabber into operation. Here is a much condensed version 
of that procedure with some helpful hints to ease the way. 
1. Getting the Board Ready 
Figure 2-1 on page 2-2 of the Frame Grabber Manual shows the jumper lo- 
cations on the board and Table 2-1 on page 2-3 tells how thev are configured from the 


factory. There are three position functions to be concerned with: 


6 Trademark of IBM 
7 Trademark of Imaging Technology Inc. 
8 Trademark of SONY 
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e Since an EGA card is needed to run “Menu” the memory base address must be 


changed on the Frame Grabber board. Address D0000 is recommended since this 
is the default value in the software. To change to address D0000 simply insert a 
jumper at location J8. If other than D0000 is chosen it should be noted for entry 
each time “Menu” 1s run. 


“Menu” has a default register base address of 100. To change the FrameGrabber 
board insert a jumper at location JI]. Again, if 100 is not used, the chosen value 
will have to be entered each time “Menu” is run. 


If an AT host is being used a jumper is required at location J20. If an XT or PC 
is being used no jumper should be installed at J20. 


2. Installing the Board 
Pages 3-1 through 3-3 of the FrameGrabber Manual provide a detailed de- 


scription for installing the board. A summary is as follows: 


Turn off and unplug the computer and periphials. 

Remove the cover. 

Place the board in any open slot where it fits and install screws. 
Reinstall cover. 


Plug in and turn on computer. 


3. Cable Installation 


Page 3-5 of the Frame Grabber Manual shows how to connect the cable. Here’s 


the bottom line: 


Plug the rectangular female connector on the cable into the top connector 
protruding from the board. 


Connect the white labeled BNC connector to the video source output. 


Connect the green labeled BNC connector to the ‘line’ input (not the ‘green’ input 
of the red, green, blue connections). 


Do not connect the red- or blue-labeled cables. 


4. Running the Diagnostics. 


Page 3-10 of the Frame Grabber Manual describes the procedure for running 


the diagnostic software included with the Frame Grabber. Since the factorv configura- 


tion has been changed as described above, the command <configure> will have to be 


entered followed by the register base address (100) and the memory base address 


(D0000). The diagnostics take several minutes to run, so be patient. As stated in the 


manual, system memory or other devices mapped into the same memorv space or regis- 


ter addresses will cause some tests to fail. 
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D. . LETS TRY IROUs 
This is an introduction to the “Menu” program. It lets you walk through some of 
the features of the program; sort of a hands-on familiarization. Turn on your source 
(camera or VCR) and video monitor now to let them warm up. On your computer go 
to the directory where you want to keep the program and copy ‘menu.exe’ there. 
I. Load the Program 
To get started simply tvpe the word ‘menu’ followed by <enter>. No <en- 
ter> is needed for any selection from here on out unless otherwise noted. The menu 
should now be displaved on the computer monitor. If vou are not using the memory 
base address (DO000) and register base address (100) as described in paragraph C above, 
hit ‘K’ to change setup and follow the instructions on your monitor. 
2. Live Action 
Menu selection ‘E’ provides continuous acquisition and display of the incoming 
real time video images. Hit ’E’ and you should see live action from your camera or VCR 
(if vou are using a tape, make sure vou press the play button first). 
3. Take a Snapshot 
Now let’s freeze the action. Hit ’F’ for snap. This acquires and displavs a single 
frame. You should see a stop action snapshot on your video monitor. The image 1s also 
stored now in framegrabber memory A, the default memory chosen when vou started 
up. 
4. Take Two Snapshots 
Now hit ‘D’ for display memory B. Take another snapshot with the ‘F’ key. 
The second snapshot is now displayed on the monitor and stored in framegrabber 
memory B. Trv going back to memory A (the ‘C’ key) to make sure the first snapshot 
is still there. Note that we are using the display memory kevs and not the select memory 
kevs to perform these operations. More on select memory later when we start storing 
images. 
5. All Black 
Let's erase one of our snapshots. Display memory B (the ‘D’ key) and then hit 
the ‘G’ key for screenclear(0). This clears the screen and framegrabber memory B to 
black or intensity level Zero. 
6. All White 
While still in memory B hit the ‘H’ key for screenclear(255). This erases the 
screen and memory B to white or intensity level 255. Darker gray shades have lower 


intensit¥ levels while lighter shades have higher intensity levels. 
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7. Put Your Snapshot in the Photo Album 

Saving an image to disk is analogous to putting vour snapshot in a photo al- 
bum. You should still have an image in frame memory A. If not, use the display mem 
A and snap routines to get one. Now hit the ‘A’ key to select memory A. This selects 
memory A for data manipulations or storage. For our purposes we have selected the 
image data in memory A to store on disk. Now hit the ‘I’ key and you will be asked to 
enter an image name. Type a name (e.g., ‘image’) to store the image on the current disk 
and directory or path and name (e.g., ‘a:image’) to store it in a different place. Then 
youll be prompted for a comment. You may give it a label or simply type ‘none’ (note 
that <enter> 1s required for each of these inputs). The disk drive selected should make 
some noise and your snapshot is safely tucked away in the “photo album” (1.e., the image 
file is written to the disk). 

8. Open the Album to See Your Snapshot 

Let’s see if the image was really saved. First erase both frame memories as 
described in paragraph 5 or 6 above. Then select a memory using ‘A’ or ‘B’. Now hit 
the ‘J’ key. Type the image file name (followed by <enter>) to read the image from 
disk. The image should appear on the video monitor. If it does not, display the memory 
you selected (the ‘C’ or ‘D’ key) and the image will appear. 

9. Simple Processing 

Now let’s process real time video images using some simple lookup table mod- 
ification techniques. Hit the “L’ key - its called analysis, but reallv it just lets you see 
how vour image is being changed. The graph vou see on your computer monitor is the 
current lookup table. The horizontal coordinate represents input pixel values and the 
vertical represents output values. 

On the top of the computer monitor screen are the current options available for 
selection. The snap (‘F’) and grab (’G’) options work the same as on the previous menu. 
Try them! Now try the invlut (‘1’) option and see how the lookup table changes (the 
blue line is the current lookup table and the green is a reference linear lookup table line). 
The effect on the image is just like looking at the negatives of a photograph, except that 
we sull have live action. Now hit the ‘L’ key for linear lookup table. We are back to 
unprocessed live images. 

Keys ‘P’ (logarithmic) and ‘A’ (absolute value) also modify the lookup table. 
You can experiment with these functions to see the effect on the real time image and 


their corresponding lookup table shape. 
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10. Histogram Equalization 

Now for the most powerful of our lookup table modification routines. Hit ’E’ 
for equalize. You now are given two choices, default and rectangle. This determines 
which portion of the image will be used to do the processing. Most of the time you 
will probably use the default mode but let’s try rectangle now just for fun. Hit ’R’ for 
rectangle. On the video monitor you will see the live images being displayed with a 
rectangle overlaid. Now hit “Num Lock’ on vour keyboard. Using the ‘2’, ‘4’, 6’, and 
‘8’ kevs on the number pad (the ones with the arrows) you can move the rectangle 
around on the screen. The corner numbers ('1’, ‘3’, ‘7’, and’9’) are used to change the 
size of the rectangle. If your image had a large white or black area in which you were 
not interested, it would lessen the impact of the equalization routine on the area of in- 
terest. So you would frame the area of interest with the rectangle. After vou have 
placed the rectangle where you want it, hit ’S’ for stop. 

Now you have to choose between manual or automatic update. Try manual 
update first. Hit the “M’ kev. The image will equalize, the new lookup table is plotted, 
and live action continues using the new lookup table. In this mode, as you notice the 
contents of the image changing you can hit the ‘U’ key and equalization will occur again. 

Let’s go back now and trv the default rectangle parameters. Hit ‘C’ to change 
parameters. Now hit ‘D’ for default and ‘C’ for continuous update. You are prompted 
for time between updates in seconds. Enter anv number 2 through 9. The computer 
will perform a new equalization at the interval vou specify. This is intended for situ- 
ations where the contents of the image are changing often and vou don’t want to have 
to remember to update. When vou have seen enough updates hit ’S’ for stop. 

11. Before and After 

Let's look at some histograms and see whats happening when we equalize. Hit 
the ‘H’ key for histogram. Again, vou may choose between default or rectangle to 
specify a portion of the image. Hit ‘D’ for default. Now hit the ’L’ key for computation 
of linear lookup table histogram. After the "COMPUTING HISTOGRAM” indicator 
goes away, hit “M’ for computation of modified histogram. When the indicator goes 
away this time, hit ‘D’ to display the linear lookup table histogram. This is the infor- 
mation used to perform the equalization. The image on the video monitor is the un- 
processed version. Now hit ’E’ to display the modified histogram. The video image is 
enhanced and the “spread out” histogram is displayed. As vou toggle back and forth 


between “D° and ‘E’. you can see graphically how equalization lets us use the full range 


of intensity levels and the corresponding enhancement of contrast on the image. Con- 
vinced? Hit the ‘S’ key and let’s trv something else. 
12. Modifying the Lookup Table 

Hit the “N’ kev. Now we'll go again to the number pad with ‘Num Lock’ on. 
Using the corner keys (‘I’, ‘3’, ‘7’, and “9°) we'll change the slope of the lookup table by 
moving the top or bottom of the curve. To shallow out the lookup table use the ‘9’ key 
to move the top to the nght and the ‘1° key to move the bottom to the left. Notice that 
enhancement to the live images changes as you modify the curve. When you are finished 
hit ‘S’ to stop and the resulting lookup table 1s shown. If you were to develop a “fa- 
vorite’ lookup table which you think you might want to use again, there is a way to save 
it in the computer. We'll talk about that later. 

13. Draw Your Own Lookup Table 

Now hit “M’ to perform another form of manual lookup table modification. 
We'll do this by drawing lines. Using the number pad with “Num Lock’ on move the 
pointer on the computer monitor by hitting the ‘2’, 4’, ‘6’, and ‘8’ keys. This will move 
it one pixel in the direction of the arrows on the keys. The next key clockwise from each 
of these will move it in the same direction ten pixels . (E.g., The ‘1’ kev moves the 
pointer in the same direction as the ‘2’ kev only faster.) When the pointer is at a lo- 
cation where you want the line to begin hit ‘B’ for begin line. Then move the pointer 
to a location where vou want the line to end and hit ’E’ for end line. Now hit ‘D’ to 
draw the line and your change will take effect on your image. You may draw other lines 
to connect segments and build a whole new lookup table. To get rid of the clutter just 
hit ‘S’ and then reenter the routine by hitting “M’. 

Now hit ‘P’ to put a pointer on the live video image. To move the pointer use 
the number pad kevs in the same manner as vou did to move the pointer on the com- 
puter monitor described above. Notice on vour computer monitor that the position of 
the pointer and the intensity level are displaved. One possible use of this capability 
would be to find the intensity level of an object in the image and use the draw lookup 
table routine to change all pixels with that intensity to another level. Again, if vou were 
to draw a reallv good all-purpose lookup table, vou could save it with a routine which 
we'll describe later. Hit ‘S’ to exit the video monitor pointer routine and ‘S’ again to exit 


the modify routine. 
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Finally, let's go back to the main menu to see the lookup table storage’selection 
routines and try some single frame (frozen) processing. Type ‘S’ to exit the analysis 
function. 

i4. Storing and Retrieving Your Favorite Lookup Tables 

Remember that lookup table you designed using the modify routine a while ago? 
Here’s how you can save it for future use. Hit the ‘M’ key. Type in a name (e.g. ‘mylut’ 
or ‘a:mylut’), hit enter, and the values are saved. To get it back hit the ‘N’ key to read 
the lookup table values from disk. Again, type in the file name followed by enter and 
vour lookup table will take effect. 

15. Piggyback Lookup Tables 

The framegrabber board has two sets of lookup tables, one as the image 1s put 
into frame memory and one as the image 1s output to the monitor. Up to this point we 
have onlv used the output lookup table. Now we'll modify both sets to achieve a dif- 
ferent result. Hit the ‘P’ key to select the input lookup table. Then go back to the 
analysis routine (‘L’) and do a histogram equalization. Everything still looks the same 
as the last ume we equalized, including the video image, but this time the computed 
values went into the input lookup table. Now exit the equalization and analvsis routines 
to get back to the main menu. Hit the ‘O’ key to select the output lookup table. Then 
go back into analysis. You may now use anv of the lookup table modification routines 
(e.g. invlut, abslut, etc.) and you will get a “piggvbacked” lookup table effect on the 1m- 
age. To get back to normal, hit ‘L’ to linearize the output lookup table. Go back to the 
main menu and select and linearize the input lookup table. Note that when vou use this 
procedure, always put values derived by equalization into the input lookup table or the 
equalization will be based on already modified data. 

16. Zoom 

Now hit the ‘R’ key to start the zoom function. On the video monitor you can 

see the rectangle we saw before. Using the number pad with ‘Num Lock’ on you can 


move the rectangle with the ‘2’, ‘4’, ‘6’, and °8’ keys in the direction of the arrows on the 
kevs. The size of this one is fixed at approximately one quarter of the screen, since we 
will be expanding the image to four times the original. Move the rectangle to the por- 
tion of the image vou wish to zoom in on and hit ‘S’. The expanded image 1s displayed. 
Notice that some resolution is lost from the original because we are “zooming” digitally 
as Opposed to changing the focal length of the camera lens. Hit “M’ to get back to the 
menu. Now you can process the “Zoomed” image using the fixed image processing rou- 


tines described below. 
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17. Unzoom 
To return to a normal image after using the zoom routine and performing any 
fixed image processing, hit the ‘R’ kev. It is essential that this be done when vou are 
finished with a “zoomed” image or the zoom parameters will interfere with other rou- 
tines. 
18. Sharpen 
Snap an image by hitting the ‘F’ key or by using the zoom routine (the ‘R’ key). 
Then hit ‘S’ to sharpen it. If you watch closely you can see the sharpening occurring 
on the displaved image. 
19. Low Pass Filter 
Snap or zoom another image and then hit ‘T’ to call the low pass filter routine. 
This function is designed to remove high frequency noise from a fixed image. Again, 
looking closelv you can see the change occurring. 
20. An Edge Detector 
One more time, snap or zoom an image. This time hit the ‘U’ key to perform 
a Sobel Edge Filter. Note the interesting display of the image’s edges on the video 
monitor. 
21. Image Averaging 
Hit the “V’ key. Then enter the number of image frames you want to average 
together. This number may be limited by the amount of disk space you have available 
(2 is a good number to start with). When the average image is displayed, vou should 
notice that objects which passed across the field of view during the operation have faded 
awav. The more frames you average, the less visible these objects will become. Image 
averaging is very effective at eliminating transient noise, but requires a stable platform 


for the source. 


To exit the program hit the “‘W’ key. 


FE. QUICK REFERENCE LIST 
This list presents the program’s options in the order they appear with a brief de- 


scription of their output. 
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1. A for select mem A 
Selects framegrabber memory A for storing or reading an image to/from disk. 
2. B for select mem B 
Selects framegrabber memory B for storing or reading an image to/from disk. 
3. C for display mem A 
Selects framegrabber memory A for acquisition of image from camera and dis- 
plays memory A on the monitor. 
4. D for display mem B 
Selects framegrabber memory B for acquisition of image from camera and dis- 
plavs memory B on the monitor. 


5. E for grab 
Continuously acquires image frames from source and displavs on video monitor. 


Uses framegrabber memory designated by ‘display mem’ function. 
6. F for snap 
Acquires a single frame from source and displays it on the video monitor. Uses 
framegrabber memorv designated bv ‘display mem’ function. 
7. G for sclear(Q) 
Clears screen and memory to intensity level 0. 
8. H for sclear(255) 
Clears screen and memory to intensity level 255. 
9. J] to store image 
Stores contents of selected memory to disk. Prompts input of filename (includ- 
ing path if desired) and comment. 
10. J to read image 
Reads image from disk file to selected framegrabber memory. Prompts input 
of filename (including path if desired). 
11. kK to change setup 
Allows change of software memory base address and register base address. On 
screen instructions and example inputs provided. 
12. L for analysis 
Displays lookup tables and histograms. Performs equalization and lookup table 
modification routines. Options as follows: 
a. m=modlut 
Modifies lookup table by drawing lines. Places pointer on video monitor. 


Move pointer to desired start point and end point using “Num Lock’ on and number pad 
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keys ‘2’, ‘4’, 6’, and °8’ for down, left, up, and right respectively. Use keys ‘1’, ‘7’, ‘9’, 
and’3’ for faster movement down, left, up, and right respectively. 

(1) b=beginline. Marks beginning of line. 

(2) e=endline. Marks end of line. 

(3) d=drawline. Draws line and enters new values into lookup table. 

(4) p=pointer. Places pointer on video monitor. Gives pointer location 
and pixel intensity level of that location on computer monitor. Use same method to 
move pointer as for line drawing routine above. 

(5) s=stop. Exit lookup table modify routine. 

b. n=mod2 lut 

Changes slope of existing lookup table and implements results. With ‘Num 
Lock’ on, use number pad keys ‘1’ and ‘3’ to move bottom of curve left and right re- 
spectively. Use ‘7’ and ‘9’ kevs to move top of curve left and right respectively. Hit ‘s’ 
to exit this routine. 

c. h=histogram 

Computes and displays histograms for images processed with linear and 
modified lookup tables. Snaps images before and after processing and allows for com- 
parison. 

(lj) d=default. Selects default portion of image for histogram computa- 
tion (approximately 90% of image). 

‘2) r=rectangle. Displays rectangle superimposed on image on video 
monitor. Change portion of image used for histogram computation by altering rectangle 
size and location. With “Num Lock’ on, use number pad kevs ‘2’, ‘4’, 6’, and 8’ to move 
rectangle down, left, right, and up respectively. Use ‘l’ and ‘3’ keys to increase and 
decrease height. Use ‘7’ and ‘9’ keys to increase and decrease width. Hit ‘s’ to exit rec- 
tangle routine. 

(3) l=com linhist. Computes histogram of image processed by linear 
LUT. Displays "COMPUTING HISTOGRAM” message until done. 

(4) m=com modhist. Computes histogram of image processed by 
modified LUT. Displays “COMPUTING HISTOGRAM” message until done. 

(5, d=dis linhist. Displays histogram computed above using image 
processed by linear lookup table. Unprocessed image shown on video monitor. 

(6) e=dis modhist. Displays histogram computed above using image 
processed by modified lookup table. Processed image shown on video monitor. 


(7) S=Stop. Exit histogram routme. 
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d. e=equalize 

Performs histogram equalization with update on command or automatically 
at specified intervals. Displays resulting lookup table after each update. 

(1) d=default. Selects default portion of image to use for equalization 
(approximately 90% of image). 

(2) r=rectangle. Allows use of rectangle to specify portion of image to 
be used as described in histogram rectangle routine above (paragraph 12.c.(2)). 

(3) c=continuous update. Selects automatic update of lookup table at 
specified intervals. Prompts for desired interval. Hit ‘s’ to exit equalization routine from 
this mode. 

(4) m=manual update. Selects lookup table update on command. 

(a) u=update— Performs new equalization and updates lookup 
table. 
(b) c=changeparam— Returns to default/rectangle portion of 
image options for changing parameters. 
(c) f=freezeframe— Acquires and displays a single frame (snap). 
(d) g=grab— Continuous acquisition and display of image 
frames. 
(ec) s=stop— Exit equalization routine from this mode. 
e. l=linlut 
Sets lookup table to linear values. Results in unprocessed image. 
fi t= invlut 
Sets lookup table to inverse values. Produces a negative image. 
g. p=loglut 
Sets lookup table to logarithmic values. Increases image brightness. 
h. a=abslut 


Inverses intensity levels greater than 127. Reduces image brightness. 


1 g=grab 
Continuous aquisition and display of image frames. 
J. f=freeze 


Aquisition and display of a single frame (snap). 
k. s=stop 


Exit analvsis routine. 


13. M to store lookup table 
Stores contents of selected lookup table to disk. Prompts input of filename 
(including path if desired). 
14. N to read lookup table 
Reads lookup table values from disk to lookup table they were stored from. 
Prompts input of filename (including path if desired). 
15. O to select output lookup table 
Selects the output lookup table to perform modifications on and for display in 
the analysis routine. When used in conjunction with ‘P’ (select input lookup table) al- 
lows for processing through two modified lookup tables. 
16. P to select input lookup table 
Selects the input lookup table to perform modifications on and for display in the 
analysis routine. When used in conjunction with ’O’ (select output lookup table) allows 
for processing through two modified lookup tables. 
17. Q for zoom 
Enlarges chosen area of image. Provides moveable rectangle to select area for 
down, left, right, and up respectively. Hit ‘Zz’ to stop rectangle and fix,zoom image. Hit 
‘m’ to return to main menu. 
18. R for unzoom 
Returns image to normal state. 
19. S for sharpen 
Sharpens fixed image using a high pass filter. May be used on zoomed image. 
20. T for low pass 
Performs a low pass filter operation on a fixed image. Reduces high frequency 
noise. Blurs the image. May be used on zoomed image. 
21. U for edge detector 
Performs Sobel edge detector operation on a fixed image. Displays detected 
edges. May be used on zoomed image. 
22. WV for image average 
Averages specified number of image frames and displays results. Prompts input 
of number of frames to average. 
23. W for framegrabber off and exit program 


Terminates program. 
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APPENDIX B. PROGRAM LISTING 


/*This program provides image processing options in a menu format. It 
was written to be compiled using Microsoft C 5.0. It also uses include 
files and the library of functions from the ITEX PCplus Image Processing 
Package. The Imaging Technology PCplus Framegrabber must be installed 
to run the program. Graphics in the program require an EGA card and 
monitor. */ 


/*List include files*/ 

include <itexpfg.h> /*from ITEX PC plus*/ 
##include <stdtyp. h> /*from ITEX PCplus*/ 
#include <stdio.h> 

#include <stdlib.h> 

##include <process. h> 

f#finclude <time. h> 

##Hinclude <conio. h> 

f##Finclude <graph. h> 


/*Declare variables*/ 

int 

i,err,n,m,t{ 10] ,s{ 10] ,j} ,highcut,r,array| 256] ,pp,qq,intens.c,a.b.da,dpe 
e,color,xl,x2,x3,x4,yyl,y2,y3,y4,mmm,nnn,res,l,low,high,blx,bly,elx,ely, 
ppl,qql,d, lowcut,start,range,intensity,y,g,h,array2[ 256] ,vt,nun; 
unsigned regbase; 

unsigned int grp; 

unsigned long membase; 

char name[ 20] ,comment[ 200] , lutname[ 20] , lutname2[ 20]; 

short grmode; 

long histvalsl{ 256]; histvals2[ 256]; histvals[ 256] ; 

FILE “stream; 

time_t ltime,xx; 

float dely,delx,f,af,bf,gf,hf,yi,delta; 

struct videoconfig config; 


/*Main program*/ 
main() 


/*Set up the framegrabber board*/ 
sethdw(0x100,0xDOO00L, DUAL); 
setdim@Si2 ,51258); 
fgon( ); 
initialize(); 
setvmask(0x00); 
sclear( 255); 


/*Initial rectangle and LUT parameters*/ 
a=30; b=30; da=450; db=450; grp=GREEN; 
loop: 


/*Present menu on screen*/ 
_clearscreen(_GCLEARSCREEN); 
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pimiatat (type: \n.); 

printf("A for select mem A\n"); 
printf(''B for select mem B\n"'); 
printf("C for display mem A\n'); 

DEintt aD tor display mem Bin"); 

PrintiG@wheeor grabin, ); 

printf("F for snap\n"); 

printf("G for sclear(0)\n"); 

printt( H for sclear(255)\n'); 

printf("I to store image\n"); 

printf("J to read imagejn'); 

printf(''K to change setup\n =) 

printf("L analysis\n'); 

printf("M for store lut\n"); 

printf("'N for read lut\n"); 

printf("0 for select output lut\n"); 

printf("P for select input lut\n"); 

printf£(''Q for zoom\n' ); 

printf( R for unzoom\n ); 

printf("S for sharpen\n'); 

printf("T for low pass\n’'); 

printf('"'U for edge detector\n"); 

printf(''V for image averagej\n' ); 

printf("W for framegrabber off and exit program\n’’); 


/*Get menu choice*/ 
i=getch(); 
i=tolower(i); 
_clearscreen(_GCLEARSCREEN); 


/*Select memory option*/ 
if(i=="a'){ 
select_mem(MEM_A); 
goto loop; } 


/*Select memory option*/ 
else if(i=='b'){ 
select_mem(MEM_B); 
goto loop; } 


/*Display memory option*/ 
else if(i='c ){ 
display_mem(MEM_A); 
goto loop; } 


/*Display memory option*/ 
else if(i=='d )§ 
display_mem(MEM_B); 
goto loop; } 


/*Continuous acquisition and display option*/ 
else if(i='e'){ 
grab(~1); 
goto loop; } 


/*Acquire and display one frame option*/ 
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else if(i=='f' ){ 
snap( WAIT); 
goto loop; } 


/*Clear screen to black option*/ 
else if(i=='g'){ 
sclear(0); 
goto loop; } 


/*Clear screen to white option*/ 
else if(i=='h'){ 
sc learG255))- 
goto loop; } 


/*Save image to disk option* / 
else Glas ie) | 
printf(’ ‘\n\n filename?\n"); 
scanf(' ‘%s" name); 
printf(’ ‘comment ?\n" ); 
scanf(''%s'', comment); 
if((err==saveim(0,0,511,511,EIGHT_BIT,name , comment ))<0){ 
switch(err) { 
case FILE_ERROR: 
printf("error opening file\n"); 
break; 
case FORMAT_ERROR: 
printf("unknown format selected\n"); 
break; 
default: 
printf("unknown error%d\n"); 
break; } 
exit(1);} 
goto loop; } 


/*Read image from disk option*/ 
else sg ioe 
prinere ‘\n\n filename?\n"); 
scanf("%s",name); 
if((err=-readim(0,0,511,511,name, comment) )<0) { 
switch(err) { 
case FILE_ERROR: 
printf("error opening file\n"); 
break; 
case FORMAT_ERROR: 
printf("unknown file format\n"); 
break; 
default: 
printf(‘‘unknown error\n"); 
break; } 
exit i je} 
printf('%s\n",comment); 
goto loop; } 


/*Change memory and register base addresses in software option*/ 


else if(i=='k )j 


regbase=0x100; membase=O0xDOO00OL; 

printf(’ ‘\n\n Default register base etdeess 460x100. \n"); 
printf("Do you want to change it? type 'y' or 'n'\n"); 
nnn=getch( ); 

if(nnn=='n _)goto nextl; 

rl y')f 

printf("Enter new register base address. e.g. 100 \n''); 
scanf("%x'' ,&regbase); } 

nextl: 

printf("\n\n Default memory base address is 0xDOOOOL. in): 
printf("Do you want to change it? type ‘y' or ‘'n'\n"); 
res=getch(); 

if(res=='n goto next2; 

if(res==' y'){ 

printf("Enter new memory base address. e.g. D0000\n"); 
scanf("%1x'',&membase); } 

next2: 

sethdw(regbase,membase,DUAL); 

goto loop; } 


/*LUT analysis _option ce | 
else if(i=='1')§ 
_setvideomode( 16); 
select_mem(MEM_A); display_mem(MEM_A); grab(-1); 
loop160: 
_clearscreen(_GCLEARSCREEN); 


/*Display current LUT*/ 
be GOoGrdt ), 
ralut(erp,0,0s256, array); 
displaylut(); 


/*Present LUT analysis options*/ 
—Settextposition(1, Ds 
printf(''m=mod lut; n=mod2; h=histogram; e=equalize; 1l=linlut; 
i=inviut une 
Settextposition(2, sp 
printf("a=abslut; p=loglut; g=grab; f=freeze; s=stop; ae 


/*Get LUT analysis choice*/ 
c=getch(); 
c=tolower(c); 


/*Continuous acquisition and display option*/ 
if(c=='g'){select_mem(MEM_A); 
display_mem(MEM_A); grab(-1); goto loop160; } 


/*Acquire and display a single frame option*/ 
else if(c=='f'){select_mem(MEM_A); 
display_mem(MEM_A); snap(WAIT); goto loop160; } 


/*Exit analysis option” / 
else if(c== s )goto stop30; 


/*Histogram option” ‘/ 
else if(c=='h )goto part2; 
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/*Histogram ,equalization option*/ 
else if(c= e ){histequal(); goto loop160; } 


/*Load linear LUT option*/ 
else if(c=='1'){linlut(grp,0); goto loop160; } 


/*Load inverse LUT option*/ 
else if(c=='i'){invlut(grp,0); goto loop160; } 


/*Load absolute value LUT option*/ 
else if(c=='a'){abslut(grp,0); goto loop160; } 


/*Load logrithmic LUT option*/ 
else if(c=='p'){loglut(grp,0); goto loop160; } 


/*Modify LUT using slope change option*/ 
else if(c=='n')goto mod; 


/*Modify LUT by drawing option*/ 
else if(c=='m')goto partl; 


else{goto loop160; } 


/*LUT drawing routine*/ 

parti: 

pp=180; qq=170; blx=pp; elx=b1x-1; 

loop1s0; 

_settextposition(1, De 

printf('b=beginline; e=endline; d=drawline; p=pointer; s=stop; 
_settextposition(2, Is 

Drantfl( 

_moveto(pp-2,qq);_lineto( pp+2 ,qq); 
_moveto(pp,qq-2);_lineto(pp,qqt2); 
c=getch(); 
c=tolower(c); 
if(c== 3") {pp=pp+10; if(pp>304 )pp=304; } 
else if(c== 7'){pp=pp-10; if( pp<51)pp=51; } 
else if(c=='1' ){qq=qq+10; if(qq>299)qq=299; } 
else if(c=='9' ){qq=qq-10; if( qq<46)qq=46; } 
else if(c= ){pp=ppt1; 4 £(pp>304)pp=304; } 
else if(csos )ipp—pp- 1, 1i(ppsow epost 
else if(c== ) {qq=qqt1; if( qq>299)qq=299; } 
else if€c= 8 ){qq=qq-1; if(qq<46)qq=46; } 
else if(c= b') {b1x=pp; bly=qq; } 
else if(c=='e' ){elx=pp; ely=qq; } 
else if(c=='s')goto loop160; 
else if(c=='d')goto draw; 
else if(c='p'){pointer(); grab(-1);} 
elsef{goto loop150; } 
goto loop150; 
draw: 
if(blx==elx)elx=bl1x+1; 
_moveto(blx,bly);_lineto(elx,ely); 
for( 1l=b1x; 1<=elx; 1=1+1) § 

dely=(float)bly-( float)ely; delx=(float)elx-( float)b1x; 


—=-_ =<—-| —_-_ —-_ =-=-.—d—ClCcrtllh lc rll 
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f=(300-( float )bly)+((dely/delx)*(1-(float)b1x)); 
array[ 1-50 =(int)f; } 

walut(grp,0,0,256,array); displaylut(); 

goto partl; 


/*Slope changing routine*/ 
mod: 
-Settextposition(1, Ip): 
printf("s=stop 
_Settextposition(2, 18 

Eammen( 

1=0; 

modi: 

if(array{ 1+1] >0){a=1l; goto mod2; }1=1+1; goto modl; 

mod2: 

1=255; 

mod3: 

mee eyl ts 1] <255) {b=1; goto mod4; }1=1-1; goto mod3; 

mod 

a g=b; for(1=0; 1<=255; 1=14+1)array2/[ 1] =array{ 1]; 

mod5: 

c=getch(); 

c=tolower(c); 

if(c=='7' )E=s- 1s 

else if(c== 9! ')g=gt1; 

else if(c=='1' )h=h-1; 

else if(c=='3' )h=h+1; 

else if(c=='s')goto loop160; 

else{goto loop160; } 

hf=(float)h; af=(float)a; bf=(float)b; gf=(float)g; 

yi=(255*(hf-af))/(bf-af-gfthf); 

for( y=0; y<h; y=y+1) array/[ y] =0; 

for(y=h; y<(int )yi; y=y+1) { 
delta=((hf-af)/(Cyi-hf))*(yi-(Cfloat)y); 
array[ y] =array2[(Cint)((float)y-delta)]; } 

EOuCy (allt) yt yoy yays 1) { 
dedeasGeet “bf Cyi-gf))*(yi-( pie jy); 
array[ y] =array2[ (int)((float)y-delta)]; } 

for( y=g; y<=255; y=y+1)array[ y] =255; 

displaylut();walut(grp,0,0,256,array); goto mod5; 


/*Histogram routine*/ 
paree. 
_settextposition(1,1); 
printf("d=default; r=rectangle; 
Searles toosition(2, ye 
ove sony Cm 
c=getch(); 
c=tolower(c); 
if(c=='d' ){ a=40; b=45; da=410; db=405; } 
else if(c=='r'){mrectangle(); grab(-1);} 
else{goto part2; } 
linlut( INPUT,0); linlut(grp,0); 
select_mem(MEM_A); display_mem(MEM_A); snap(-1); 
walut( INPUT,0,0,256, array); 
select_mem(MEM_B); display_mem(MEM_B); snap(-1); 


By 
Ne 


ye 
>: 


linlutC INPUT ,0); 
next45: 
_settextposition(1,1); 
printf("l=com linhist; m=com modhist; 
d=dis linhist; e=dis modhist;"); 
_Ssettextposition(2, 1); 
printf(''s=stop a 
c=getch(); 
c=tolower(c); 
if(c=='1' )comvalsi1(); 
else if(c='m' )comvals2(); 
else if(c=='d'){_clearscreen(_GCLEARSCREEN); lutcoord( ); displaylut(); 
histcoord(); plotvals1(); linlut( GREEN ,0); display_mem(MEM_A); } 
else if(c=='e'){_clearscreen(_GCLEARSCREEN); lutcoord(); 
displaylut(); histcoord( ); plotvals2(); walut(GREEN,0,0,256,array); 
display_ mem(MEM_ Wie 
else if(c=='s'){linlut(GREEN, 0); walut(grp, 0,0,256,array); 
_settextposition(2,1);print£(" ny. 
select_mem(MEM_A); display_mem(MEM_A); grab(- lye 
goto loop1é0; } 
else{goto next45; } 
goto next45; 


/*Exit analysis*/ 


stop30: 

Settextposition(1, itp 

printeecs Ee 
_Settextposition(2, ily is 

printf( Bo 


system(''mode C080"); 
a=30; b=30; da=450; db=450; 
goto loop; } 


/*Save LUT to disk option*/ 
else ine m'){ 
Princes enter LUT file name\n”); 
scanf("'%s'"', lutname); 
save_lut(lutname,grp,0,0,256,NEW_FILE); 
goto loop; } 


/*Read LUT from disk option*/ 
else as een 
printf(” enter LUT file name\n"); 
scanf(''%s"', lutname2); 
read_lut(lutname2); 
goto loop; } 


/*Select output LUT for loading option*/ 
else if(i='o'){ 
grp=GREEN; 
Force Loop ..} 


/*Select input LUT for loading option*/ 


else if (i=='p'){ 
grp=INPUT; 
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goto loop; } 


/*Zoom option*/ 
else if(i=='q' ){ 
select _mem( MEM_ A); display_mem(MEM_A); 
print £C ‘\n\n To move upper left corner of rectangle type: \n\n"); 
printf("' 8" for up\n"); 
Drintt@s2 for down\n ); 
primes + for leftin ); 
printf("'6' for right\n\n\n"); 
printf("type 'z' to zoom\n\n'); 
loop/0: 
while(! kbhit()){snap( WAIT); rectangle(a,b,240,235,0); } 
e=getch(); 
e=tolower(e); 
if(e=='8') {b=b-10; if( b<0)b=479; } 
else if(e=='2'){b=b+10; if(b>479)b=0; } 
else if(e='4" ) {a=a-10; i f(a<0)a°5125 } 
else if(e= ee _) {a=atl0; if(a>512)a=0; } 
else if(e=='z')goto stop3; 
elsef{goto loop/7Q; } 
goto loop/0; 


ttt 


stop3: 

Sey eo MICL), zoom(HAG.. Zo 

spot20: 

printf("\n\n type 'm' to return to menu\n’); 
ee 


c=tolower(c); 
if(c=='m' ){da=240; db=235; goto loop; } 
else{goto spotz20; } 


/*Unzoom option*/ 

else if(i='r'){ 
zoom(MAG_1);scroll(-b); pan(-a); 
select_mem(MEM_B); display_mem(MEM_B); sclear( 255); 
select_mem(MEM_A); display_mem(MEM_A); sclear( 255); 
initregs();setlut(grp,0); 

select_mem(MEM_A); display_mem(MEM_A); 

gmap <1); 

a=40; b=45; da=410; db=405; 

Soto Loop; } 


/*Sharpen option*/ 


else Le Claas Si 
printf("\n\njn\n\ n\n\n\n\n\n\n 
PROCESSING IMAGE\n"); 
sharpen(a,b,da,db, POSITIVE); 
goto loop; } 


/*Low pass filter option*/ 
else if(i='t')f{ 
printf("\nin\n\nin\n\n\n\n\n\n 
PROCESSING IMAGE\n"); 
lopass(a,b,da,db, POSITIVE); 
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goto loop; } 


/*Sobel edge filter option*/ 


else if(i=='u'){ 
printf(‘'\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING IMAGE\n''); 
sobel(a,b,da,db,16); 

goto loop; } 


/*Image average option*/ 


else if(i='v')f{ 
printf("number of frames?\n'); 
scanf(''%d" ,&num); 
_clearscreen(_GCLEARSCREEN); 
printf("\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING IMAGE\n"); 
average(10,10,490,470,num,num); 
goto loop; } 


/*Framegrabber off option*/ 


else iis) 
CHORDIE De 


/*This function places a moveable pointer on video screen and feeds 
position and intensity level to computer monitor*/ 
pointer(){ 


select_mem(MEM_A); display_mem(MEM_A); 

_settextposition(1,1); 

printf(''s=stop age 
_settextposition(2,1); 

prune cee ae) 
ppl1=255; qq1=240; 

loop50: 

snap(WAIT); 

intens=rpixel({ppl,qql); 

_settextposition(2,65); 

printf(' x=%d;n'',pp1); 

_settextposition(3,65); 

printf("y=%din" ,qq1); 

_settextposition(4,65); 

printf("intensity=%d\n' ,intens); 

if(intens<=127)color=255; 

ifCintens?127)color=0; 

x1=pp1-4; x2=pp1+4; yyl=qq1; y2=qq1; x3=pp1; x4=pp1; y3=qq1-4; y4=qq1+4; 
clr_line(xl,yyl,x2,y2,color); cir_line(x3,y3,x4557-, core). 
while(! kbhit()){snap(WAIT); clr_line(xl,yy1,x2,y2,color); 
clr_line(x3,y3,x4,y4,color); } 

c=getch(); 

c=tolower(c); 

if(c==' 3° ){ppl=pp1+10; if(pp1>511)pp1=0; } 

else if(c== 7'){ppl=pp1-10; if( pp1<0)pp1=511; } 

else if(c=='1'){qql=qq1+10; if(qq1>479)qqi=0; } 

else if(c=='9'){qql=qq1-10; if(qq1<0)qq1=479; } 
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else 
else 
else 
else 
else 
else 


if(c=='6'){ppl=pp1+1; if(pp1>511)pp1=0; } 
if(c== Sys See 1; if(pp1<0)pp1=511; } 
if(c=='2') {qql=qqltl; if(qq1>479)qq1=0; } 
if(c=='8'){qql=qql-1; if(qq1<0)qq1=479; } 
if(c=='5' )goto loop50: 


if(c=='s')goto stop; 


elsef{goto loop50; } 


goto 


loop50; 


StOp:. 
—Settextposition(2, may 
prints 
—Settextposition(3, cs), 
printt( 
_Settextposition(4, 65)," 
primer. ea co 


changed on the video screen. 
used in other functions*/ 
mrectangle()§ 


/*This function places a rectangle whose size and location may be 
Size and location parameters are then 


select_mem(MEM_A); display_mem(MEM_A); 
_Settextposition(1, Wie 

printf(''s=stop 

_settextposition(2,1); 

Print t( 

loop60: 


while(! kbhit()){snap(WAIT); rectangle(a,b,da,db,0); } 


e=getch(); 
e=tolower(e); 


if(e=='8'){b=b-10; if (b<0)b=479; } 

else if(e=='2' ){b=b+10; if(b>479)b=0; } 

else 1 amabaey mbes aa ig 

else if(e=='6' ){a=at10; if(a>512)a=0; } 

else if(e=='9' ){da=dat10; if(data>512)da=512-a; } 
else if(e=='7' ){da=da-10; if(da<10)da=10; } 

else if(e=='3' ){db=db+10; if(db+b>479)db=479-b; } 
else if(e=='1' ){db=db-10; if(db<10)db=10; } 

else if(e=='s')goto stop2; 

else{ goto loop60; } 

goto loop60; 

Spool: 


histequal(){ 


loop3: 

_settextposition(1,1); 
printf('‘d=default; r=rectangle; 
_settextposition(2,1); 

prance( 

c=getch(); 

c=tolower(c); 

if(c=='d' ){a=40; b=45; da=410; db=405; } 
else if(c=='r'){mrectangle(); grab(-1); } 
else{goto loop3; } 
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/*This function performs histogram equalization and allows for manual or 
automatic updates*/ 


"ys 
oe 


leopinl: 
_settextposition(1, iy 
printf(‘'c=continuous update; m=manual update; 
_settextposition(2, Lye 
print re: 

c=getch(); 

c=tolower(c); 

if(c=='c ‘goto loop12; 

else if(c=='m'); 

elsef{goto loop1]; } 

color=4; 

loop2: 

select_mem(MEM_B); display_mem(MEM_B); 
snap(WAIT); 

select_mem(MEM_A); 
display_mem( MEM_A); 

grab(-1); 

select_mem( MEM_B); 
histeq(grp,0,a,b,da,db); 
ralut(grp,0,0,256,array); 

“setcolor(color ye tlie 
displaylut(); 

color=color+t5; if(color>9)color=4; 

loop4: 

_settextposition(1,1); 
printf("u=update; c=change param; f=freezeframe; g=grab; 
_settextposition(2,1); 
piaAnerGe 

c=getch(); 

c=tolower(c); 

if(c=='u')goto loop2; 

else sale c'){_setcolor(9 % 16); goto loop3; } 
else if(c== ae _d{snap(WAIT); goto loop; } 

else if(c== 4S "){grab(-1); goto loopd4; } 

else if(c=='s'){_setcolor(9 % 16); goto stop25; } 
else{goto loop4; } 

loop12: 

color=4; 
_settextposition(1,1); 
printf('enter time between updates (integer 2-9) 
_settextposition(2,1); 
Pritcte 

c=getch(), 

if(c=='2' )vt= 2; QS shee" o )vt=3; else if(c=='4' )vt=4; 
else if(c=='5' \vtss; else if(c=='6')vt=6; 


Ne 
ap i: 


s=stop;"’); 


woe 


me) 
eG 


else if(c=='7')vt=7;else if(c=='8' )vt=8; else if(c=='9' )vt=9; 


else{goto loop12; } 
_settextposition(1,1); 
printf("c=change param; s=stop; 


_settextposition(2,1); 
primcLes 
while(! kbhit() ){ 
time(&ltime); 
xx=ltime; 
select_mem(MEM_B); display_mem(MEM_B); 


50 


OE 


a 


snap(WAIT); 
select_mem(MEM_A); 
display_mem(MEM_A); 
grab(-1); 
select_mem(MEM_B); 
histeq(grp,0,a,b,da,db); 
Palitieemp. 00 ,256,array ); 
_setcolor(color % 16); 
displaylut(); 
color=colort5; if(color>9)color=4; 
shide¢(itime-xx)<vt)time(&ltime); } 
c=getch(); 
c=tolower(c); 
if(c='c'){_ setcolor(9 % 16); goto loop3; } 
else if(c='s'){_setcolor(9 % 16); goto stop25; } 
else{_setcolor(9 % 16); goto loop3; } 
stop25: 
select_mem(MEM_A); 
display_mem(MEM_A); 
} 


/*This function displays the current LUT shape on the computer monitor*/ 
displaylut()§{ 

c=_getcolor(); 

_setcolor(2 % 16); 

_moveto(50,300);_lineto( 305,45); 

_setcolor(c % 16); 

_moveto(50,300); 

for( 1=50; 1<=304; 1=1+1) { 

_lineto(1+1,300-array[ 1-49] );3 3 


/*This function draws the coordinates for the histogram display*/ 
histcoord() { 
_moveto(370,300);_lineto(370,40);_moveto(370,300);_lineto( 630,300); 
_settextposition( 23,47); print £("%d", oe 
Bocrueweseoition(2), 56). printf(' %d ,100): 
_settextposition( 23,71); printf('%d", 200); 
Beettextosition( 15,43), printf(''%d',2); 
_settextposition(8,43); printf("%d", (i ))R 
for( 1=370; 1<=620; 1=1+50) { 
“moveto( 1,300); _lineto(1,305); 
_moveto(370,670-1);_lineto(365,670-1);} } 


/*This function computes a histogram for the contents of memory A*/ 
comvals1(){ 
_Settextposition(1, 1p 
prineee ao) 
_settextposition(2,1);printf(" : 
_settextposition(13,50); printf(' ‘COMPUTING HISTOGRAM’); 
select_mem(MEM_A); 
histogram(a,b,da,db,1,1,0,histvalsl); 
_settextposition( 13,50); printf(" oD) sat 


/*This function computes a histogram for the contents of memory B*/ 
comvals2() { 
_settextposition(1,1); 


prone ia =): 
_settextposition( 2,1); printf(" So 
_settextposition( 13,50); printf(''COMPUTING HISTOGRAM"); 
select_mem(MEM_B); 

histogram(a,b,da,db,1,1,0,histvals2); 

_settextposition(13,50); printf(" ye 


/*This function plots the comvalsl output histogram*/ 
plotvals1(){ 
_moveto( 370,300); 
for( 1=370; 1<=625; 1=1+1) § 
_moveto( 1,300); 
_lineto(1,300-(histvals1[ 1-370] /20)); }3 


/*This function plots the comvals2 output histogram*/ 
plotvals2(){ 
_moveto( 370,300); 
for( 1=370; 1<=625; 1=1+1){ 
_moveto( 1,300); 
_lineto(1,300-(histvals2[ 1-370] /20)); }} 


/*This function draws the coordinates for the LUT display*/ 
lutcoord( ) { 
_moveto(50, 300); __lineto(50,40); _moveto( 50,300); _lineto(310, 300); 
_settextposition( 23,7); printf('%d",0); 
_settextposition( 23,18); printf(''%d'', 100); 
_settextposition( 23,31); printf("%d"', 200); 
Uscttextposition( 15,3): pranta@ecc OOD. 
_settextposition(8,3); printf('"%d" , 200); 
for( 1=50; 1<=300; 1=1+50) § 
_moveto( 1,300); _lineto( 1,305); 
_moveto(50,350-1);_lineto(45,350-1); }} 
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